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Nomenclature 


BILCE  Bivariate  Integration  LCE 
Ck  Coefficient  of  order  K  Hermite  polynomial 
CLT  Central  Limit  Theorem 
CR  Capability  ratio 
DOL  Day  of  launch 

fx  Density  function  of  a  random  variable  x 

Hk  Order  K  Hermite  polynomial 

Lai  low  Structural  allowable  load,  lb  or  in-lb 

Laxiai.  Axial  load,  lb  or  in-lb 

Lgy  Equivalent  axial  load,  lb  or  in-lb 

LCE  Loads  combination  equation 

Mp  Pitch  bending  moment,  lb  or  in-lb 

Mt  Combined  yaw  and  pitch  bending  moment,  lb  or  in-lb 

My  Yaw  bending  moment,  lb  or  in-lb 

m,  Third  central  moment 

N(|l,a)  Normal  random  variable  with  mean  p  and  standard  deviation  0 
P  Random  variable  representing  Mp 
Pr  Probability 
RSS  Root  sum  of  squares 

U(a,b)  Uniform  random  number  with  minimum  a  and  maximum  b 
X  Random  variable  representing  Mx 

Y  Random  variable  representing  My 

Px  Mean  of  random  variable  X 

Ox  Standard  deviation  of  random  variable  X 


5 


6 


1.  Introduction 


Most  launch  vehicles  can  only  achieve  the  desired  level  of  structural  reliability  by  restricting  the 
winds  through  which  the  vehicle  is  allowed  to  fly.  This  restriction  is  accomplished  by  first 
analytically  flying  the  vehicle  through  wind  profiles  that  are  measured  just  prior  to  launch,  and 
calculating  altitude  histories  of  angles  of  attack,  dynamic  pressure,  rigid  body  acceleration,  and 
engine  gimbal  angles,  among  others.  These  time  histories  are  then  used  to  establish  static-aeroelastic 
and  other  day-of-launch  calculated  loads,  which  are  then  combined  with  the  pre-day-of-launch 
calculated  loads  to  obtain  the  total  load.'  "  Since  many  of  the  individual  loads  have  random 
characteristics,  a  statistical  load  enclosure  of  the  total  load  is  computed.  This  total  load  enclosure  is 
then  compared  to  the  vehicle  allowable  strength,  and  if  it  is  exceeded,  the  vehicle  is  not  launched.  If 
sufficient  time  is  available  before  the  end  of  the  launch  window,  the  entire  process  is  repeated  for 
subsequent  wind  profiles.  If  there  is  not  enough  time,  the  launch  attempt  is  aborted,  and  the  vehicle  is 
prepared  for  the  next  available  launch  window. 

For  the  launch  vehicle  core  structure,  it  is  typical  to  express  the  comparison  between  the  day-of- 
launch  calculated  load  enclosures  and  the  allowable  values  as  capability  ratios  (CRs);  i.e., 

_  Equivalent  Axial  Load  _  ( 1 ) 

StructuralAllowableLoad 


The  structural  allowable  load,  L^llow.  is  determined  prior  to  the  day  of  launch,  usually  with  structural 
testing  and  appropriate  factors  of  safety.  The  equivalent  axial  load,  L^y,  is  determined  from  an 
evaluation  of  the  stress  at  a  cross  section  of  the  launch  vehicle.  It  is  a  function  of  the  resultant 
bending  moment,  M^,  obtained  from  the  combination  of  the  yaw  and  pitch  bending  moment 
components.  My  and  Mp.  It  also  depends  on  the  axial  load,  L^xial-  ^nd  the  radius  to  the  point  on  the 
launch  vehicle  cylinder  wall; 


L  =L 

^EQ.mX  ^AXIALMAX  ^  „  J- 

Radius 


(2) 


The  combined  pitch  and  yaw  bending  moments,  Mj,  are  computed  by  combining  a  set  of  statistically 
varying  loads  from  a  variety  of  sources.  To  ensure,  with  a  high  probability,  that  the  predicted 
equivalent  axial  load  is  not  exceeded  during  atmospheric  flight,  a  load  enclosure  of  Mj  is  used  in  Eq. 
2  to  compute  the  CRs.  The  load  enclosure  is  usually  computed  at  a  99.7  percent  enclosure,  90 
percent  confidence,  level. ^ 

The  enclosure  load  is  computed  by  combining  loads  obtained  from  analyses  performed  to  simulate 
the  various  atmospheric  flight  load  contributors.  Static-aeroelastic  (STEL)  loads  analyses  are 
performed  to  establish  the  loads  that  are  due  to  that  portion  of  the  vehicle’s  angle  of  attack  that  vary 
relatively  slowly  with  time.  The  STEL  load  is  a  function  of  the  day-of-launch  winds  and  the  vehicle 
steering  profile.  Gust  analyses  are  performed  to  establish  the  response  of  the  vehicle,  and  its  payload, 
to  the  turbulence  that  might  be  encountered  on  any  given  flight.  Buffet  analyses  are  performed  to 
establish  the  loads  due  to  the  dynamic  response  of  the  vehicle/payload  system  to  shock  waves,  flow 
separation  due  to  changes  in  vehicle  geometry,  and  the  interaction  between  the  two.  In  addition, 
other  analyses  are  often  performed  to  estimate  loads  due  to  items  such  as  vehicle  load  alleviation 
steering,  autopilot  noise,  wind  measurement  error,  changes  in  day-of-launch  winds  from  the  time  they 
are  measured  to  when  the  vehicle  is  to  be  launched,  and  vehicle  dispersions  from  the  nominal 
parameters  used  in  the  analyses. 

The  most  accurate  method  of  computing  the  enclosure  load  is  by  Monte  Carlo  simulation.**  Due  to 
time  constraints,  Monte  Carlo  simulation  is  currently  not  a  realistic  computational  option  during  day- 
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of-launch  operations.  An  analysis  and  launch  decision  is  required  within  a  matter  of  a  few  minutes, 
and  a  Monte  Carlo  simulation  requires  considerably  more  time.  In  order  to  quickly  calculate  CRs  on 
the  day  of  launch,  a  computationally  efficient  analytical  procedure  is  necessary. 

A  widely  used  analytical  method  for  computing  the  enclosure  load  uses  what  are  referred  to  as  Loads 
Combination  Equations  (LCEs).'^'^  Enclosure  loads  and  the  resulting  CRs  computed  using  the  LCEs 
have  been  shown  to  be  biased  on  the  conservative  side.  This  bias  reduces  launch  availability. 
Reduction  factors  have  been  incorporated  into  the  LCEs  to  reduce  the  conservative  bias.  However, 
this  can  also,  on  occasion,  lead  to  underpredicting  the  CRs. 

This  paper  presents  the  derivation  of  a  new  procedure,  referred  to  as  the  Bivariate  Integration  Load 
Combination  Equation  (BILCE),  that  has  less  bias  and  variance  than  other  LCE  procedures. 
Comparisons  of  the  performance  of  the  BILCE  and  LCE  are  made  by  using  both  to  derive  load  CRs 
from  several  heavy  lift  launch  vehicle  load  datasets,  and  comparing  these  to  load  CRs  computed  using 
a  Monte  Carlo  simulation. 
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2.  DOL  Implementation  of  LCE 


An  example  of  an  actual  LCE  input  load  dataset  from  the  day  of  launch  of  a  heavy  lift  launch  vehicle 
is  given  in  Table  I .  Only  the  STEL,  lack-of-wind  persistence,  and  wind  measurement  error  loads  are 
established  on  the  day  of  launch,  since  these  depend  on  the  day-of-Iaunch  winds  and  steering  profile. 
The  remaining  loads  are  computed  prior  to  the  day  of  launch.  Most  of  the  loads  are  assumed  to 
follow  a  normal  distribution,  except  the  buffet  and  autopilot  loads,  which  are  assumed  to  follow  a 
Rayleigh  distribution,^  and  the  gust  load,  which  is  assumed  to  follow  a  gamma  distribution.’ 


Table  1.  Predicted  Bending  Moments  (Million  Inch-Pounds) 


Yaw  Plane 

Pitch  Plane 

Source 

Computation 

Distribution 

Mean 

Standard 

Mean 

Standard 

Deviation 

Deviation 

STEL 

DOL 

- 

97.53 

0 

635.61 

0 

Wind  Persistence 

DOL 

Normal 

86.11 

141.87 

68.95 

109.37 

Wind 

DOL 

Normal 

0 

82.59 

0 

91.08 

Measurement 

Error 

Trajectory 

Pre-DOL 

Normal 

0 

100.18 

0 

96.28 

Dispersion 

Aerodynamic 

Pre-DOL 

Normal 

0 

3.16 

0 

19.91 

Maneuvering 

Pre-DOL 

Normal 

167.49 

24.56 

147.55 

19.86 

Wind  Data  Gap 

Pre-DOL 

Normal 

0 

0 

0 

0 

Root  Sum  of  Squares 

Root  Sum  of  Squares 

Buffet 

Pre-DOL 

Rayleigh 

119.72 

141.22 

Auto  Pilot 

Pre-DOL 

Rayleigh 

0 

0 

Response  to  30-fL/sec  Gust 

Response  to  30-ft/sec  Gust 

Gust 

Pre-DOL 

Gamma 

1976.58 

1909.37 

A  pictorial  of  how  the  enclosure  load  is  determined  is  given  in  Figure  1 .  In  the  figure,  IMy  and  EMp 
are  the  yaw  and  pitch  vector  sums  (summed  means  +  RSS  of  dispersed  terms),  while  R  is  the  99.7 
percent  enclosure  load.  The  steps  to  compute  the  enclosure  load  are:  1)  Compute  the  mean  and 
dispersed  term  of  each  load  in  the  yaw  and  pitch  directions.  The  dispersed  term  is  the  load’s  99.7 
percentile  minus  its  mean.  2)  Reduce  the  means  and  dispersed  terms  of  the  gust,  buffet,  autopilot, 
trajectory,  and  maneuvering  loads  by  ten  percent,  as  a  typical  value.  This  is  the  reduction  factor 
obtained  from  a  Monte  Carlo  simulation  (Ref.  1).  3)  Add  the  sum  of  the  means  (excluding  STEL) 
with  the  RSS  of  the  dispersed  terms  in  the  yaw  and  pitch  planes  (Fig.  la).  4)  Translate  the  yaw  and 
pitch  vector  sums  to  the  STEL’s  yaw  and  pitch  vector  load  (Fig.  lb).  5)  Construct  an  ellipse  such  that 
major  and  minor  axis  vertices  are  the  translated  yaw  and  pitch  vector  sums  (Fig.  Ic).  6)  The 
enclosure  load  is  the  magnitude  of  the  point  on  the  ellipse  that  is  the  maximum  distance  from  the 
origin  (Fig.  Id). 
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Figure  1.  Illustrates  the  LCE  method  to  compute  the  limit  load.  ZMy  and  ZMp  are  the  yaw  and 
pitch  vector  sums  (i.e.,  summed  mean  +  RS  of  dispersed  terms).  R  is  the  99.7  percent 
enclosure  load. 


Further  details  of  the  methodology  used  in  this  LCE  can  be  found  in  Ref.  I .  The  LCE  has  sufficed 
for  many  years  and  many  launches  by  various  heavy  and  medium  lift  launch  vehicles.  Although  it 
has  served  its  purpose,  the  approach  has  a  certain  inherent  bias  and  variability  that  are  due  to  the 
following  reasons;  1 )  The  mean  sum  added  to  the  RSS  of  the  dispersed  terms  (steps  3  and  4)  is  not  a 
99.7  percent  enclosure  for  the  sum  of  non-normal  loads.  An  example  to  illustrate  this  is  the  sum  of 
two  independent  U(0,1)  random  variables.  The  99.7  percentile  of  the  sum  is  1.92,  whereas  the  mean 
sum  added  to  the  RSS  of  the  dispersed  terms  is  1.70.  2)  The  magnitude  of  the  point  on  the  ellipse 
that  is  the  maximum  distance  from  the  origin  (steps  5  and  6)  is  not  the  99.7  percent  bound  of  the 
load  magnitude,  even  if  the  vertices  of  the  ellipse  are  the  true  99.7  percent  load  bounds.  An  example 
here  is  the  RSS  of  two  independent  N(0,1)  random  variables  where  the  center  of  the  ellipse  is  the 
origin.  The  99.7  percentile  of  the  RSS  is  3.41,  whereas  the  magnitude  of  the  line  from  the  origin  to 
the  point  on  the  ellipse  furthest  from  the  origin  is  3.89.  Even  with  these  approximations,  Monte 
Carlo  simulations  have  shown  that  the  LCE  generally  performs  as  intended.  However,  it  is  now 
possible  to  combine  the  various  load  contributors  in  a  more  rigorous  manner.  It  is  the  purpose  of 
this  paper  to  introduce  such  a  procedure. 
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3.  DOL  Implementation  of  the  BILGE 


The  enclosure  load  of  the  bending  moment  is  the  value  R  such  that 

<  /?)  =  0.997  (3) 


where  Y  and  P  are  random  variables  representing  the  sum  of  the  yaw  and  pitch  bending  moments,  My 
and  Mp.  Given  that  Y  and  P  have  density  functions,  fy  and  fp,  Eq.  (3)  can  be  rewritten  as 

JJ  fYMfp(v)dudv  =  0.997  (4) 

Vh’ +!■-<« 


Due  to  the  non-normal  loads  (buffet,  autopilot,  and  gust)  and  the  numeric  difficulties  in  convoluting 
their  density  functions  with  the  other  bending  moment  densities,  computing  fy  and  fp  and 
consequently  the  exact  value  of  R  in  a  timely  manner  during  day-of-launch  operations  is  time 
consuming. 

Since  Y  is  the  sum  of  several  statistically  independent  loads,  from  the  Central  Limit  Theorem,  its 
density  function,  fy,  is  approximately  normal  with  mean,  py,  the  sum  of  the  individual  load’s  means, 
and  standard  deviation,  Oy,  the  RSS  of  the  individual  load’s  standard  deviations.  This  approximation 
is  sufficient  as  long  as  the  dominating  loads  are  approximately  normal.  Unfortunately,  two  significant 
load  contributors,  buffet  and  gust,  have  Rayleigh  and  gamma  distributions,  respectively,  which  are 
both  highly  skewed  to  the  right.  The  Central  Limit  Theorem  approximation  can  be  improved"  by 
modifying  the  normal  density  approximation  so  that  its  higher  central  moments  are  equal  to  the 
higher  central  moments  of  Y.  The  modified  density  approximation  can  be  written  as 


exp 


i  +  lQW* 

i  =  l 

(5) 


where  H^fx)  is  a  Hermite  polynomial  of  order  k,  and  Ck  are  coefficients  determined  by  equating  the 
central  moments  of  Y  with  the  density  approximation.  Equating  only  the  first  three  moments,  we  find 


C,  =  C,  =  0  and  C,  =  -4 
'  ^  6al 


(6) 


where  m,  is  the  third  central  moment  of  Y;  i.e.,  [E(Y-py  )'].  Since 


Hyix)  =  x^-3x 


(7) 


Eq.  (5)  can  be  written  as 
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1 


(8) 


fyM^ 


j,V2^ 


•  1  + 


3_ 

6ol 


orj 


The  approximate  density  function  of  P  is  determined  in  an  analogous  manner.  Using  ^s.  (4)  and 
(8),  the  steps  to  compute  the  BILGE  enclosure  load  are:  1)  Compute  the  mean,  standard  deviation, 
and  3"*  moment  (skewness)  of  each  load  in  the  yaw  and  pitch  direction.  These  statistics  m^  be 
computed  analytically  or  by  Monte  Carlo  techniques  for  the  pre-day-of-launch  loads.  2)  Sum  the 
means,  variances,  and  3"*  moments  to  determine  p.,  a,  and  m,.  3)  Approximate  the  density  functions, 
fv  and  fp,  of  the  bending  moment  sums,  Y  and  P,  using  Eq.  (8).  4)  The  enclosure  load  is  the  value  R 
in  Eq.  (4)  found  by  using  Romberg  numerical  integration  (Ref.  2)  and  the  approximate  density 
functions,  fy  and  fp. 

Though  higher  moments  can  be  used  to  approximate  fy  and  fp,  for  the  launch  vehicle  flights 
examined,  the  first  three  moments  were  adequate  to  account  for  the  skewness  as  a  result  of  the  buffet 
and  gust  loads. 
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4.  Accuracy  of  the  BILCE 


The  accuracies  of  the  LCE  and  the  BILCE  are  assessed  by  comparing  their  CRs  with  CRs  computed 
using  a  Monte  Carlo  simulation.  Details  of  the  Monte  Carlo  simulation  are  given  in  Refs.  3,  4,  and  9 
Comparisons  are  made  at  several  launch  vehicle  stations  and  flight  Mach  numbers  from  nine  heavy 
lift  launch  vehicle  flights. 

A  comparison  of  the  CRs  for  one  of  the  flights  at  Mach  1 .58  is  shown  in  Fig.  2.  The  marks  that 
overlay  the  top  horizontal  line  are  CR  differences  (Monte  Carlo  CR  -  LCE  CR)  derived  at  48  launch 
vehicle  stations.  The  bottom  marks  are  CR  differences  at  the  same  stations  using  the  BILCE. 


Figure  2.  Marks  indicate  the  difference  between  the  Monte  Carlo  CR  and  the  LCE  and 
BILCE  CRS  at  48  launch  vehicle  stations  at  mach  1.58  for  a  heavy  lift  launch 
vehicle.  The  improved  accuracy  of  the  BILCE  compared  to  the  LCE  can  be  seen. 

The  reduction  in  variation,  using  the  BILCE  compared  to  the  LCE,  demonstrates  the  improved 
accuracy  of  the  method.  For  the  BILCE,  the  lack  of  points  far  to  the  left  of  zero  indicates  fewer 
overpredictions  of  the  enclosure  load  compared  to  the  LCE,  providing  the  potential  for  making  fewer 
unnecessary  “No-Go”  launch  decisions  by  using  the  BILCE.  The  lack  of  points  far  to  the  right 
indicates  the  potential  to  also  make  fewer  “Go”  launch  decisions  with  lower  than  the  desired 
statistical  enclosure. 

The  CR  comparisons  for  the  same  flight  at  the  other  Mach  numbers  studied  are  given  in  Fig.  3. 
Results  are  similar  to  the  Mach  number  1.58  values;  CR  accuracy  is  improved  using  the  BILCE.  CR 
comparisons  for  the  remaining  eight  heavy  lift  launch  vehicle  flights  show  similar  results. 
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Figure  3.  Additional  CR  difference  plots  at  various  mach  numbers  (Fig.  2  included  on  right  in  second  row 
from  top).  Marks  indicate  the  difference  between  the  Monte  Carlo  CR  and  the  LCE  and  BILCE 
CRS  at  48  launch  vehicle  stations  for  a  heavy  lift  launch  vehicle.  Each  plot  shows  the  improved 
accuracy  the  BILCE  compared  to  the  LCE. 

14 


5.  Conclusions 


A  refined  procedure  for  combining  atmospheric  flight  loads  to  compute  load-to-allowable  ratios  has 
been  presented.  It  is  shown  that  this  procedure  produces  CRs  with  less  bias  and  variation  than  a 
previously-used  procedure.  As  a  result,  day-of-launch  implementation  of  the  proposed  procedure 
should  increase  launch  availability  without  reducing  launch  reliability. 
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